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A test was conducted in the 15 cm x 15 cm supersonic wind tunnel at NASA Glenn Research 
Center that focused on corner effects of an oblique shock-wave/boundary-layer interaction. 
In an attempt to control the interaction in the corner region, eight corner fillet configurations 
were tested. Three parameters were considered for the fillet configurations: the radius, the 
fillet length, and the taper length from the square corner to the fillet radius. Fillets effectively 
reduced the boundary-layer thickness in the corner; however, there was an associated penalty 
in the form of increased boundary-layer thickness at the tunnel centerline. Larger fillet radii 
caused greater reductions in boundary-layer thickness along the corner bisector. To a lesser, 
but measureable, extent, shorter fillet lengths resulted in thinner corner boundary layers. 
Overall, of the configurations tested, the largest radius resulted in the best combination of 
control in the corner, evidenced by a reduction in boundary-layer thickness, coupled with 
minimal impacts at the tunnel centerline. 


Nomenclature 

a = speed of sound, m/s 

Hi = incompressible shape factor 

/ = fillet length, cm 

M = Mach number 

p = static pressure, kPa 

Pu = static pressure at the farthest upstream centerline static tap, kPa 

r = fillet radius, cm 

sa = uncertainty in speed of sound, m/s 

sM = uncertainty in Mach number 

sTs = uncertainty in static temperature, K 

su = uncertainty in velocity, m/s 

t = fillet taper length, cm 

Ts = static temperature, K 

u = velocity in the streamwise direction, m/s 

5 = boundary-layer thickness, cm 

5c = boundary-layer edge location measured along the corner bisector from the right-angle corner, cm 

I. Introduction 

W HEN a shock wave interacts with the low-momentum fluid in a viscous boundary layer, the pressure rise 
associated with the shock wave causes the boundary layer to thicken or, for stronger shocks, causes flow 
separation. This interaction occurs in a wide range of high-speed flows such as flow over transonic airfoils and flow 
in propulsion flowpaths. 

Because the presence of these interactions is common, much attention has been devoted both to understanding 
the physics of such flows as well as developing means to control the effects of the interaction. This research generally 
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takes the form of simplified unit problems where the phenomena can be studied in a more easily understood 
environment. Some common unit problems that have been used to investigate shock-wave/boundary-layer 
interactions include the swept fm problem^ impinging oblique shock interactions^’"^, compression ramp 
interactions^’^, normal shock interactions^’^, and a symmetric wedge compression corner^’^^. 

Many of these unit problem experiments are conducted in small wind tunnels with rectangular cross sections. 
There are several reasons for using a rectangular cross section, including ease of model integration and change out, 
ease of incorporating optical diagnostic techniques, and similarity to certain flight geometries. However, when a 
rectangular cross section is used, it introduces the complication of comers where two viscous boundary layers interact. 

A number of subsonic^ and supersonic experimental and computational studies have focused on 
understanding the flow field in an axial right angle corner, and have shown that corner flows are inherently three- 
dimensional. In the absence of transverse pressure gradients, corner secondary flow of Prandtl’s second kind 
develops^^, with velocities vectoring toward the wall along the corner bisector and away from the corner along the 
walls, as shown in Figure 1. 



Figure 1. Schematic of secondary flow structure in the corner. 

In the present study, an oblique shock wave was generated and reflected off of the floor of a small scale supersonic 
wind tunnel with a square cross section, as shown in Figure 2. This configuration incorporates three distinct unit 
problems: the swept fm interaction on the sidewalls, an impinging/reflecting oblique shock- wave interaction on the 
floor, and the comer flow interaction. 



Figure 2. Three-dimensional representation of an inviscid 

APPROXIMATION OF THE REFLECTING OBLIQUE SHOCK INTERACTION. 


Looking first at the interaction on the sidewall of the tunnel, the swept fm interaction causes the flow to turn 
parallel to the shock generator incidence angle and drives additional low momentum fluid from the sidewall boundary 
layer toward the floor of the wind tunnel and into the corner region. Because of the subsonic layer on the sidewall, 
this flow turning begins upstream of the shock wave. 

A second interaction along the floor of the tunnel is the impinging shock that spans the tunnel width and subjects 
the flow to the pressure rise of the shock and its reflection. In a viscous flow, the pressure rise occurs over some 
interaction length. 
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The interaction of the shock- wave with the comer flow is the focus of this study. Along the corner bisector, the 
subsonic layer is large compared to the boundary layers along the centerline of the tunnel span. The swept fin 
interaction introduces an additional mass of low momentum fluid to the corner region, as noted earlier. Additionally, 
as with the impinging shock interaction, the comer flow region is subjected to a strong streamwise adverse pressure 
gradient over a short streamwise extent. 

These three effects of all three interactions - the additional low moment fluid from the sidewall boundary layer, 
the large pressure rise of the shock and reflection, and the already large subsonic layer in the corner - combine to 
cause the flow in the corner to separate earlier than either the side or bottom walls resulting in a highly three- 
dimensional flow field^^. This recirculation extends somewhat upstream of the centerline interaction region, and the 
flow experiences the separation region as an aerodynamic blockage which creates a weak conic compression in the 
flow^"^’^^. The aerodynamic blockage and the compression wave affect the shape and magnitude of the impinging 
shock interaction across the entire span of the tunnel^. 

Previous studies have considered several ways to mitigate the separation in corner shock- wave interactions. Bleed 
has been successfully used to control the corner interaction^^'^^. Removing a portion of the low momentum flow from 
the corner region improves the resistance to separation. However, in one case the required flow rate of mass removal 
was too high to be considered favorable^^, and in another, the bleed showed little effect 

Aerodynamic fences have been used in supersonic wind tunnel testing where the interest is in making 
measurements on the centerline. The fences, placed parallel to the sidewalls, isolate the center region of the test 
section from the corner interaction at the tunnel sidewalls. A new, thin boundary layer develops along the fences, 
which approximates an inviscid walP^'^^. 

Vortex generators have been studied in an attempt to use mixing to increase the near wall velocity in the corner 
region to reduce corner separation, but generally have not been found to be effective at controlling the flow in 
corners^^’^^. 

Recently, chamfers applied to the corners of the flow domain in a computational study were considered as a way 
to limit the extent of the corner effects for an oblique shock interaction by blocking the migration of flow from the 
sidewall swept fm interaction into the corner. The effectiveness was shown to be similar to that in a computed bleed 
comparison case for controlling the corner interaction but avoided the need for mass removal^^. 

Fillets have also been studied as a possible means to control corner effects in many configurations with varying 
effectiveness. An early supersonic inlet study by Valerino^'^ found that incorporating fillets into the corners of the 
inlet could increase the critical total pressure recovery or extend the stable operating range depending on the Mach 
number. Additionally, tests of a subsonic serpentine inlet showed that a small radius in the corners gave the lowest 
distortions^. On the other hand, a different supersonic inlet configuration tested at low speeds found corner fillets to 
have little effects^. A computational study of a supersonic rectangular nozzle found that comer rounding of 10% of 
the nozzle half width effectively controlled the corner vortex and instability Tests of fillets on a wing body junction 
at low speed with the fillet wrapped around the entire base of the wing did not prevent separation, but merely acted to 
increase the effective radius of the leading edge of the wingS^. 

In this investigation, fillets were introduced into the comers of the wind tunnel. It was expected that the fillets 
would reduce the thickness of the subsonic layer in the corner, thereby improving the boundary-layer properties in the 
comer and increasing the resistance of the comer to separation. The geometric properties of the fillets were varied to 
determine the effect of the fillet radius, fillet length, and the length of the taper from a square corner to the fillet radius 
on the effectiveness of the fillets. 


II. Experimental Setup 


A. Facility 

The experimental testing for this study was completed in the 1 5 cm x 1 5 cm supersonic wind tunnel, shown in 
Figure 3, in the Engine Research Building of the NASA Glenn Research Center. This tunnel is a continuous flow 
facility with air supplied at 377 kPa and exhausted to a vacuum system that maintains an absolute pressure of 13.8 
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kPa. A valve upstream in the inflow line regulates the incoming total pressure. During this test, the pressure was 
controlled to maintain a Reynolds number of 13. 1E6 /m. For this test, the facility fixed wall Mach 2 nozzle was used, 
which provided a measured Mach number at the nozzle exit during the test of 1.98. 

The test section is a 1 5 cm square constant area duct extending 7 1 cm. A photograph of the test section is shown 
in Figure 4 with the x-, y-, and z-directions indicated. The coordinate system origin for this test was chosen such that 
X = 0 cm at the test section entrance, y = 0 cm on the bottom wall, and z = 0 cm on the left wall when looking 
downstream. 



Figure 3. Photograph of the 15 cm x 15 cm supersonic wind tunnel. 

The test section includes four removable wall sections, one on each wall, which can be modified as needed for 
testing. For this test, sidewall inserts were used which had large plexiglass windows. The top wall insert had a smaller 
glass window. Two different inserts were used for the bottom wall during different phases of testing. One was a 
blank, and the other was instrumented with static taps, which will be described in more detail later. 

The aft end of the test section mates with a duct that has a 23.5 cm square cross section. Incorporated in this 
section is a translation mechanism that allows for remotely controlled translation of a probe in the y- and z-directions, 
and manual translation in the x-direction. 

B. Test Hardware 

The test hardware for this study included an oblique shock generator plate and several fillet configurations. The 
shock generator and a fillet configuration can be seen installed in the test section in Figure 4. 



Figure 4. Photograph of the test hardware installed in the test section. 
(Axes shown to indicate direction, not origin location.) 
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A schematic of the layout of the test hardware in the test section is shown in Figure 5. Two representative fillet 
configurations are shown in red and green, representing the largest and smallest configurations tested. The blue dashed 
line indicates the undisturbed boundary-layer thickness measured at an x-location of 40.2 cm. The oblique shock 
generator plate, as well as the shock and the leading edge of the expansion fan it creates are shown in dotted black 
lines in the figure. The rectangle shown represents the extent of the sidewall plexiglass windows. 
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Figure 5. Schematic showing the position of the shock generator and fillet configurations. 
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Shock Generator 

To create an oblique shock-wave interaction, a flat plate with a sharp leading edge was used that spanned the 
entire tunnel. On the sides of the plate, pins extended into the plexiglass windows to allow the plate to pivot. The aft 
end of the plate was connected by a bar to an electric motor mounted in the translation section that controlled the plate 
angle. 

The location of the shock generator plate can be seen in Figure 5. The test data were collected with the shock 
generator at 8.5 degrees. However, the tunnel would not start with the plate at that angle, so during startup, the angle 
was set at approximately 6 degrees and then remotely actuated to the test setting. At 8.5 degrees, the shock generator 
leading edge location was at x = 40.0 cm, which led to an inviscid shock impingement location of x = 53.8 cm. 

The length of the shock generator plate was chosen such that it was long enough that the expansion fan from the 
aft end of the plate didn’t interfere with the measurements, but short enough so that blockage from the plate at 8.5 
degrees didn’t introduce a normal shock in the test section. 

Fillets 

The fillets for this test were created using a moldable material, EXAMIX™ NDS. A study conducted by 
Braafladt^^ looked at ten potential materials for use in this application. EXAMIX^m NDS is marketed for dental 
applications and is primarily used as an impression material. Because of this, it is designed to hold its shape within a 
tight tolerance as it cures. That property combined with a short cure time, ease of cleanup, and low surface roughness 
were the primary factors that lead to its selection. 

Initially, conducting the test with metal fillets was considered, but for the smallest radius there was concern that 
there would not be enough material for a fastener to thread into. Additionally, there was more potential for damage if 
a metal fillet were to detach during the test compared to a moldable fillet. 

The fillets were installed with a linear taper from the square corner of the test section to the full radius of the fillet. 
The radius then remained constant from that point to the end of the test section. Fillets were only applied to the two 
bottom comers of the test section as shown in Figure 6, which is where measurements were taken. 
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Figure 6. Photograph of fillets installed in the corners of the test section, looking 

DOWNSTREAM. 

Figure 7 is a close up of the fillet region from Figure 5, which shows the three geometric variables of the fillets 
that were changed during the test: the fillet radius, r; the fillet length, /; and the taper length, t. 

The fillet radii tested ranged from 0.64 to 1.59 cm. Based on boundary-layer thickness measurements at the 
tunnel centerline obtained in previous tests, these radii covered a range of approximately r!5= 1/3 to r!5= 4/3. The 
area change for the largest fillets resulted in a test section cross section that was 0.48% smaller than the square 
(baseline) cross section. This area change resulted in the test section Mach number being reduced from 1.98 to 1.97 
when fillets with the largest radius were installed. 

The fillet length, /, was measured from the start of the full radius to the inviscid shock impingement location, as 
shown in Figure 7. The fillet lengths tested ranged from 0 cm to 27.4 cm. In terms of expected boundary-layer 
thickness, this equates to a range from 1/5= 0 to 1/5= 24. 

The taper length, t, was the distance over which the corner transitioned from the square corner of the test section 
to the full fillet radius. Taper lengths of 3.4 cm and 13.7 cm were tested, representing an expected range of t!5= 3 
and t!S= 12. 


shock 




Figure 7. Illustration of the three fillet parameters that were varied. 

The fillets were formed in two sections: the taper, and the constant radius section. The constant radius section 
was formed by applying the material to the corner and using a metal cylinder of the desired radius as a mold. Tape 
was placed on the tunnel wall a distance of the radius away from the corner prior to applying the EXAMIX^m NDS. 
This allowed any excess material to be peeled away after it cured. 

The taper section was formed in the same way except that a rapid prototype mold for the transition was used 
instead of a cylinder. After both sections were formed, sand paper was used to smooth all of the edges to ensure a 
smooth transition between the two sections, as well as between the fillet and the tunnel wall. 
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C. Instrumentation 

During testing, pressure and temperature data were collected for the flow field, as well as position measurements 
and limited schlieren and oil flow visualization images. The angle of the oblique shock generator was measured with 
a rotary encoder which attached to one of the pivot pins that extended through the plexiglass window. 

The tunnel conditions during the test were recorded by instrumentation standard to the facility. The total pressure 
and temperature were measured in the plenum. The facility mass flow rate was measured with a venturi in the system 
upstream of the plenum. Two static pressures were averaged for the facility Mach number calculation. The static 
pressure taps were located on either side of the converging-diverging nozzle, near the nozzle exit. 

The bulk of the data recorded and analyzed for this effort was pressure data. Static pressures were recorded on 
the tunnel floor with 73 static taps on the centerline and 10 static taps located 2.5 cm off of the centerline on either 
side, as shown in Figure 8. 
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Figure 8. Location of static pressure taps. Each dot represents a tap location. 

A Pitot probe was mounted on the translation mechanism so that an array of total pressures could be measured 
downstream of the shock reflection. The probe was positioned with the probe tip at x = 55.4 cm, between the shock 
reflection and the leading edge of the expansion fan from the shock generator. The primary measurement plane 
location can be seen as the purple line in Figure 5. The goose-necked probe was angled at 45 degrees to be able to 
measure into the corner. The translating probe was designed with an electric circuit that closed when the probe was 
in contact with the wall. The probe position was calibrated before each run using the voltage change to determine the 
probe touch location. Where possible, a wind-on tare was completed as well. 

The distributions of points measured for the no fillet case and for the largest fillet case are shown in Figure 9. 
The same number of points was used for the fillet configurations, but the spacing was compressed. The 
measurements included 12 boundary-layer profiles: one in the corner along the bisector, seven z-profiles, and four y- 
profiles. The comer boundary-layer profile consisted of 50 points. Each of the other profiles consisted of 45 points. 
The profiles were clustered with more measurements near the wall, and data points farther apart as the points 
approached the core flow. Additionally, 27 points scattered between the profiles were collected to smooth out 
contour plots created from the data. 

D. Experimental Uncertainty 

Experimental measurements have some bias and precision errors associated with them. While we cannot know 
the exact magnitude of these errors, estimates can be made of the uncertainty associated with the measurements. 

The 103 kPa pressure transducers used for this study have a built in uncertainty of 0.05% of full scale. The as- 
recorded uncertainty, which accounts for each component in the system, increases to 0.07% of full scale, or ±0.07 kPa 
equivalently. The type E thermocouples used have an as-recorded uncertainty of ±1.4 K. 

Based on a propagation of uncertainty analysis. 

Table 1 shows calculated values and the associated uncertainty for Mach number, static temperature, the speed of 
sound, and velocity for three flow conditions. These three conditions represent a typical near wall subsonic case, the 
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supersonic boundary- layer edge case downstream of the reflected shock, and the core flow upstream of the shock from 
top to bottom respectively. 




(B) Largest fillet 

Figure 9. Array of measured Pitot pressure locations. 


Table 1. Propagation of uncertainty for three flow conditions. 
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E. Test Cases 

A total of 10 cases were included in the final test matrix for this effort. These included an upstream measurement 
plane to observe the inflow, a baseline measurement downstream of the oblique shock with no fillet, and eight fillet 
configurations. 

For the upstream profile, the probe tip was moved to x = 40.2 cm. The measurement plane location can be seen 
in Figure 5. Measurements were made at each of the locations shown for the no fillet case in Figure 9. Additionally, 
boundary-layer surveys were recorded for y = 4.5 cm, y = 6.0 cm, and y = 7.25 cm. The range of the probe did not 
reach y = 1.5 cm. These surveys mapped one quadrant of the tunnel inflow. 

The downstream baseline was measured at the same x- station as the fillet configurations. Eight fillet 
configurations were tested while varying three geometric parameters. Four levels were tested for the fillet radius and 
length, and two levels were tested for the taper length. The combinations tested are listed in Table 2. Values in the 
table have been normalized by the boundary-layer thickness measured upstream, S= 1.24 cm. 

The levels of the variables were selected to attempt to bound the problems within the confines of the available 
space in the test section. The smallest value of each variable was expected to have little effect on the flow. The largest 
value of the fillet radius was based on the lip of the schlieren frame that was available to support the fillet on the 
sidewall during formation. The largest values of length and taper were a combination which represented the farthest 
upstream location that could be accessed through the sidewall openings. 


Table 2. fillet configurations. 
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Results 


A. Upstream Flow Field 

A pressure survey was taken at a plane upstream of the interaction region in order to understand the inflow 
conditions and confirm that the inflow boundary layer was turbulent. A measurement location of x = 40.2 cm was 
selected. Additional profiles were measured along the tunnel sidewall to survey one quadrant of the flow field. 

The Mach number contours resulting from this upstream survey are shown in Figure 10, while Figure 1 1 shows 
lines of isovelocity normalized by the core velocity with the values labeled along the lines. The Mach number contours 
show a freestream Mach number near 2.0 as expected. Calculated edge Mach numbers range from 1.91 to 2.00. A 
subsonic region can be seen in the corner with lower Mach numbers than are visible elsewhere along the wall. The 
subsonic layer measured for the comer profile was nearly twice as thick as that for any other profile. This region can 
also be observed in the isovelocity contours. The lowest measured normalized velocity away from the comer was 
0.5 1 . The lowest measured normalized velocity for the corner profile was 0.35. 

Additionally, both plots show a flow pattern typical for rectangular nozzles with contouring on two walls. On the 
floor of the tunnel where the nozzle was contoured, the boundary layer was relatively uniform across the span. On 
the sidewall where the nozzle was flat, a bowing of the contours can be seen with the boundary layer thinner near the 
corner and thicker near the center of the span. 


9 



^ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ^ 

2 4 6 

z (cm) 

Figure 10. Mach contours for the upstream plane measured at x = 40.2 cm. 
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Figure 11. Iso velocity contours for the upstream plane measured at x = 40.2 cm. 

The boundary-layer thickness for each profile along the sidewall and the floor is shown in Figure 12. The abscissa 
of the plot shows positions on both the sidewall (left) and floor (right). The negative values are the locations of the 
y-profiles on the sidewall, whereas the positive values are the locations of the z-profiles on the floor. Zero is marked 
as a dashed line to indicate the location of the corner. The slight asymmetry about the corner bisector is likely due to 
the non-uniform expansion along the nozzle sidewalls (z=0) where viscous effects retard the expansion in the near- 
wall region. 

The outlier for the profile at z = 6 cm was due to a boundary-layer profile that had some slope to it even in the 
core region. Removing that point, the average boundary-layer thickness along the floor was 1.24 cm, which was used 
to normalize distances through the remainder of these results. 
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Figure 12. Boundary-layer thickness at each profile location for the upstream case (x = 40.2 cm). 

Dotted line indicates mean value along the floor. 


A similar plot can be seen in Figure 13 showing the incompressible shape factor for each profile. The shape 
factor was quite uniform throughout the flow field, and the measured values averaging 1.37 are consistent with a fully 
turbulent flow field. The comer profile had a higher shape factor, with a value of 1.58. 
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Figure 13. Incompressible shape factor at each profile location for the upstream case. 
Dotted line indicates mean value along the floor. 


B. Baseline 

A survey was completed for a baseline configuration at the primary measurement plane, x = 55.4 cm, to 
understand the effect of the oblique shock interaction on the comer flow without fillets. Photographs of oil flow 
visualization and schlieren are shown in Figure 14. The predicted inviscid shock impingement location is shown in 
the figure as a solid white line. The extent of the interaction region is shown by the dashed white lines. Multiple lines 
are shown for the interaction extent because the range is different in the comer compared to the centerline (z = 1.5 
cm). The window frame interferes with the field of view for the schlieren near the wall. The x-axis lies along the 
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tunnel wall location, so the lines tracking the shock wave angle intersect to show the location and extent of the 
interaction region along the tunnel floor. 

It can be seen that the interaction region was upstream of the calculated inviscid shock location. This was due to 
viscous effects on the oblique shock generator plate as well as the tunnel floor. The formation of a boundary layer on 
the oblique shock generator plate effectively increased the shock generator incidence angle, leading to a slightly 
steeper shock angle. The impinging interaction on the tunnel floor created a separation region that distributed the 
interaction. 



Figure 14. Oil flow visualization (top) and schlieren image (bottom) for the baseline case. 

Mach 2.0, 8.5° shock-incidence angle. 

Looking at the oil flow, the interaction region extended farther upstream in the corners. The upstream edge of 
the region across the span was straight. The downstream edge of the interaction region was curved, and the farthest 
effects downstream were seen at the centerline. Though the interaction was clearly three-dimensional, the flow field 
indicated by the oil flow was symmetric about the centerline. 

The features in the oil flow can be correlated with the features in the schlieren image. Because the schlieren 
integrates the effects across the tunnel span, this three-dimensional interaction is seen as a splitting of the shock and 
reflection in the schlieren image. The elements of the shock and reflection visible upstream correspond to the effects 
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in the corner region, and the elements visible downstream correspond to the effects in the central region of the tunnel. 
It can be seen that the extent of the interaction region agrees well between the schlieren and oil flow for both the comer 
and centerline. 

The oil flow image is repeated in Figure 15 with the shape of the interaction region indicated by an overlay. The 
region is compared to a computational solution generated by Benek et al.^ in a study looking at the effect of wind 
tunnel width on an oblique shock interaction. The computational streamlines along the floor are shown for the 
computational conditions that most closely match the present experiment. In the computational study, the width to 
height ratio was 1.0, with a thick boundary layer entering the interaction region. While the Mach number and shock 
incidence angle were different (Mach 2.9, incidence angle of 13 degrees for the computational study), the results line 
up nicely and provide a clearer picture of what may be happening in the interaction region. 



(A) Mach 2.0, 8.5° shock-incidence angle 



(b) Mach 2.9, 13° shock-incidence angle 

Figure 15. Comparison of oil flow visualization to computational results for a similar 

CONFIGURATION BY BENEK"*. 

The extent of the pressure rise in the static pressure profile along the tunnel centerline, shown in Figure 16, 
agrees well with the interaction region from the oil flow and schlieren. Ten percent of the total rise in static pressure 
had occurred by x = 46.7 cm. Ninety percent of the total rise in static pressure had occurred by x = 52.3 cm. 
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Figure 16. Static pressure distribution on the tunnel centerline for the baseline 

CONFIGURATION. 

The Mach number contours and isovelocity lines resulting from the pressure survey are shown in Figure 17 and 
Figure 18 respectively. The average edge Mach number downstream of the shock was Mach 1.44, compared to a 
calculated Mach number of 1 .38 downstream of a reflected shock at Mach 2.0 with 8. 5 -degree incidence angle. While 
the velocity was lowest in the corner region, at the measurement plane no separation was observed. 



z(cm) 


Figure 17. Mach number contours for the baseline case. 



Figure 18. Isovelocity contours for the baseline case. 


The boundary-layer thickness and incompressible shape factor for each profile are shown in Figure 19 and Figure 
20. The shock interaction increased the average boundary-layer thickness along the tunnel floor to 1.42 cm and the 
mean incompressible shape factor to 1.65. In the corner the boundary-layer thickness increased from 1.61 cm before 
the shock to 2.33 cm downstream, and the incompressible shape factor increased from 1.58 to 1.89. 
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Figure 19. Boundary-layer thickness at each profile location for the baseline case. 
Dotted line indicates mean value along the floor. 
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Figure 20. Incompressible shape factor at each profile location for the baseline case. 
Dotted line indicates mean value along the floor. 


C. Fillet Test Matrix 

Pitot pressure surveys for each of the eight fillet configurations were analyzed to calculate the velocity field 
downstream of the oblique shock at x = 55.4 cm. The lines of isovelocity for each of the cases are shown in Figure 
21 arranged from smallest fillet radius to largest. 

Compared to the baseline case from Figure 18 it can be seen that along the corner bisector, the isovelocity lines 
are spaced more closely for the fillet configurations, indicating a steeper velocity gradient near the wall. The lines 
also generally slant upward along the tunnel floor from the corner toward the centerline. For the smaller two radii, 
the shorter fillet length, /, from each pair has a steeper slope, suggesting a shorter length resulted in a more distorted 
boundary layer at the tunnel floor centerline. For the two larger radii, less difference is seen at the centerline as the 
length was changed. 
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(C) rl 6 = 0.77, H 6 = 7.35, tl 6 = 2.75 



Figure 21. Normalized isovelocity contours for the fillet configurations. 
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z (cm) 

(H) r!S= 1.28, HS= 22.0, t!S= 11.0 

Figure 21. Normalized isovelocity contours for the fillet configurations, (concluded) 

Compressing the data from the contour plots into boundary-layer thickness at each profile location, the effects of 
the fillet configurations can also be seen in Figure 22. A similar trend can be seen across the cases. Starting on the 
sidewall the boundary-layer thickness decreased approaching the corner. Because the boundary- layer profiles 
measured at 0.75 cm and 1 .5 cm on each wall were within the boundary layer of the other wall, those data points were 
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not included on the plot. Moving along the tunnel floor toward the centerline, the boundary-layer thickness gradually 
increased. 

In some cases, such as r!5= 1.28, 1/5= 22.0, t!5= 11.0 seen in part H of Figure 22, there was a peak in the 
boundary-layer thickness for the profile measured at z = 2.25 cm. Looking at the corresponding contours in part H of 
Figure 21, this can be seen as an associated bulging of the isovels in this region, indicative of a counter-clockwise 
rotating streamwise vortex. 
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Figure 22. Boundary-layer thickness at each profile location for each of the fillet 

CONFIGURATIONS. 
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Figure 22. Boundary-layer thickness at each profile location for each of the fillet 

CONFIGURATIONS. (CONCLUDED) 


Looking at equivalent plots for the incompressible shape factor at each profile location for each fillet 
configuration, shown in Figure 23, some trends can be observed. Typically, the lowest measured incompressible 
shape factor was in the corner region. For the smallest three fillet radii, the shape factor increased moving toward the 
tunnel centerline, with larger increases seen for the shorter fillet lengths. For the largest radius the levels were more 
uniform, though the largest incompressible shape factor for these cases was consistently measured at the sidewall 
profile at y = 3 cm. 
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Figure 23. Incompressible shape factor at each profile location for each of the fillet 

CONFIGURATIONS. (CONCLUDED) 


In order to compare across the cases, the boundary-layer thickness along the corner bisector is shown in Figure 
24 for all eight fillet configurations as well as the baseline. All of the fillet configurations resulted in a thinner 
boundary layer compared to the baseline. Larger fillet radii resulted in thinner boundary layers, as might be expected. 
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Looking at the effect of fillet length, for each pair of data points at a given radius {rIS) in Figure 24 the open 
symbol represents the shorter length, /. Interestingly, the shorter cases caused a greater reduction in boundary-layer 
thickness along the bisector than the longer configurations. The two cases at r!5= 1.02 have a similar magnitude. 
This may be attributed to the fact that the sum of US and tIS for the two cases was very similar, with combined lengths 
of 17.45 and 18.35 for the longer and shorter US, respectively. This was not the case for the other fillet configurations. 

Prior to the test, it was anticipated that the configurations with 1/5= 0 would have the least effect on the interaction. 
However, the benefit of these cases may in part be attributable to the local boundary layer thinning due to stream tube 
contraction resulting from the compression of the taper acting in the interaction region. This effect is minimal for the 
r/5= 0.51, 1/5= 0, t/5= 11.0 case because the small radius and long taper formed at most an angle of 1.1 degrees. 
This small angle, even using a two-dimensional approximation that is generous, would result in a local stream tube 
area reduction of 3.2 %. For the r!S= 1.28, HS= 0, t!S= 2.75 case, however, the compression angle of the taper was 
10.9 degrees, resulting in a maximum stream tube area reduction of 24 %. Because of the location of the taper, this 
effect was applied in the interaction region for the HS= 0 cases. This approximation is optimistic, but suggests that 
the effect may have locally thinned the boundary layer sufficiently to improve the flow through the interaction region. 

Examining the change in the flow field due to the fillets another way. Figure 25 shows the corner boundary-layer 
edge location versus fillet radius. The boundary layer edge location is calculated as the distance from the 90 degree 
corner to the boundary layer edge along the corner bisector. For the baseline case where the wall location is equivalent 
to the 90-degree corner, the boundary layer edge location is equal to the boundary layer thickness, S. Examining the 
results in this way shows that the physical location of the boundary layer edge changes relatively little throughout all 
of the configurations, so the thinning of the boundary layer is due to the lowest momentum region in the corner being 
filled in by the fillet. While the changes are small, it is interesting to note that for six of the eight cases, the boundary 
layer edge location has actually decreased with the introduction of the corner fillets. The largest increase in boundary 
layer edge location was measured for the r!S= 1.28, HS= 22.0, t!S= 11.0 fillet configuration, with an increase from 
2.33 cm for the baseline configuration to 2.38 cm. 

In trade for the benefits seen in the comer. Figure 26 shows the boundary layer thickening that occurred on the 
tunnel centerline (z = 7.5 cm) compared to the baseline case for all eight fillet configurations. The three smallest fillet 
radii caused similar levels of thickening at the centerline. For the smallest two radii, the length also played a factor, 
with larger centerline boundary-layer thicknesses occurring for shorter fillet lengths. Again, the magnitude of the 
boundary-layer thickness was similar for both cases with rl5= 1.02 where the combined length including the fillet 
and taper were nearly equal. On the centerline, the largest radius had the least negative impact. The boundary-layer 
thickness increased from 1.51 cm for the baseline to 1.66 for the large radius cases. Overall, of the configurations 
tested, the large radius seems to give the best combination of control in the corner, evidenced by reduction in boundary- 
layer thickness, coupled with minimal impacts at the tunnel centerline. 
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Figure 24. Boundary-layer thickness along the corner bisector versus fillet radius. 
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Figure 25. Boundary layer edge location measured along the corner bisector versus fillet 

RADIUS. 
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Figure 26. Boundary-layer thickness at the tunnel centerline versus fillet radius. 
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Looking at the incompressible shape factor for the corner profile versus fillet radius in Figure 27, the fillets caused 
a small reduction. This indicates an increased resistance to flow separation. Figure 28 shows the effect on 
incompressible shape factor at the tunnel centerline. The three smallest fillet radii produce a large increase in the 
incompressible shape factor at the centerline. Again, the largest fillet configuration shows less adverse impact at the 
centerline. In this case the longer fillets at the largest radii have the least impact with an increase in incompressible 
shape factor from 1.62 for the baseline to 1.89. 
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Figure 27. Incompressible shape factor for the corner profile versus fillet radius. 
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Figure 28. Incompressible shape factor at the tunnel centerline versus fillet radius 


Oil flow visualization was only captured for two of the fillet configurations. Figure 29 shows the resulting images 
compared to the oil flow for the baseline case. The solid line in each figure indicates the calculated inviscid shock 
impingement location. The dashed lines show the extent of the interaction region in the corners and along the 
centerline for the baseline case to allow comparison between the images. Because the oil did not flow well on the 
EXAMIX material, it is somewhat difficult to observe the effect of the fillet on the corner region. 

For the r!S= 0.51, HS= 0, t!S= 1 1.0 case, it can be seen that the upstream extent of the interaction region across 
the span has moved downstream. Additionally, the angle of the oil flow which follows the corner separation appears 
shallower. For the r!S= 1.02, HS= 14.7, t!S= 2.75 case, the oil flow patterns are very similar to those observed for 
the baseline case. 
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(B) rl6= 1.02, H6= 14.7, t!6= 2.75 



(c) Baseline 


Figure 29. oil flow visualization for two fillet configurations compared to the baseline. The 

SOLID LINE INDICATES THE INVISCID SHOCK IMPINGEMENT LOCATION. THE DASHED LINES INDICATE THE 
EXTENT OF THE SEPARATION REGION FOR THE BASELINE CASE. 
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Conclusions 

The shock- wave/boundary-layer interaction tested included three distinct unit problems: a swept fin interaction, 
an impinging oblique shock-wave interaction, and a comer interaction. The corner interaction was the focus of this 
study. 

Each of the elements of the configuration produced an effect in the corner region. The swept fm interaction turned 
the flow on the sidewall toward the floor, driving additional low momentum fluid into the corner region. The 
impinging shock interaction had a greater static pressure rise than the swept fm interaction because the oblique shock 
and its reflection encountered the floor boundary layer in the same region. The effect of the viscous wall was to 
distribute this interaction region. This pressure rise across both the shock and reflection was also experienced in the 
comer. Because the wall shear in the corner must be zero (at y = 0, du! dz = 0), the corner was in a state of incipient 
separation. This also resulted in a thick subsonic layer in the corner through which the pressure rise could propagate. 

Because of this combination of influences, the corner region was the first to separate, and the large subsonic 
region allowed the effect to extend upstream. In order to attempt to control the interaction in the corner region, corner 
fillets were applied with various radii, fillet lengths, and taper length from square to the full fillet radius. 

Fillets were observed to be effective in reducing the boundary-layer thickness in the corner downstream of an 
oblique shock. There was an associated penalty in the form of increased boundary-layer thickness at the tunnel 
centerline. The fillets effected control on the comer boundary layer through a reduction in the thickness of the subsonic 
layer in the corner, as well as by removing the enforced region of localized zero wall shear stress. Additionally, for 
the l/S= 0 cases, a benefit may have been gained from the effect of a localized reduction in boundary-layer thickness 
due to stream tube contraction from the compression created by the taper. 

Comparing the results for the eight cases, it was found that a larger radius caused a greater reduction in boundary- 
layer thickness along the corner bisector. To a lesser, but measureable, extent the shorter fillet lengths resulted in 
thinner comer boundary layers. At the centerline, boundary layer thickening was seen for all of the cases. The largest 
fillet radius resulted in the least negative impact at the centerline. Little effect of taper was noted. Overall, of the 
configurations tested the large radius seemed to give the best combination of control in the corner, evidenced by 
reduction in boundary-layer thickness, coupled with minimal impacts at the tunnel centerline. 
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